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Abstract 



We review some properties of periodic orbit families in polygonal billiards 
and discuss in particular a sum rule that they obey. In addition, we provide 
I algorithms to determine periodic orbit families and present numerical results 

o. 

\& , that shed new light on the proliferation law and its variation with the genus of 

ON 

£^ I the invariant surface. Finally, we deal with correlations in the length spectrum 

>v 

and find that long orbits display Poisson fluctuations. 
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X 1 I. INTRODUCTION 

Billiards are an interesting and well studied class of hamiltonian systems that display 
a wide variety of dynamical behaviour depending on the shape of the boundary. The ray 
equations commonly considered follow from a short wavelength expansion of the Schrodinger 
equation with Dirichlet (or Neumann) boundary conditions. A particle thus moves freely 
between collisions at the boundary where it suffers specular reflection. As an example, a 
square billiard generates regular dynamics that is restricted to a torus in phase space due to 
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the existence of two well behaved constants of motion. In contrast, generic trajectories in the 
enclosure formed by three intersecting discs explore the entire constant energy surface and 
hence the system is ergodic. Moreover, orbits that are nearby initially move exponentially 
apart with time and hence the system is said to be metrically chaotic. 

Polygonal billiards are a class of systems whose properties are not as well known. The 
ones accessible to computations are rational angled polygons and they generically belong 
to the class of systems referred to as pseudo-integrable |§. They possess two constants 
of motion like their integrable counterparts @] but the invariant surface T is topologically 
equivalent to a sphere with multiple holes || and not a torus. Such systems are characterised 
by the genus g where 2g is the numbers of cuts required in T to produce a singly connected 
region (alternately, g is the number of holes in T). As an example, consider the 1-step 
billiard in Fig. 1. For any trajectory, p\ and p 2 y are conserved. The invariant surface consists 
of four sheets (copies) corresponding to the four possible momenta, (±p x ,±p y ) that it can 
have and the edges of these sheets can be identified such that the resulting surface has the 
topology of a double torus. 

The classical dynamics no longer has the simplicity of an integrable system where a 
transformation to action and angle co-ordinates enables one to solve the global evolution 
equations on a torus. On the other hand, the dynamics is non-chaotic with the only inter- 
esting feature occurring at the singular vertex with internal angle 37r/2. Here, families of 
parallel rays split and traverse different paths, a fact that limits the extent of periodic orbit 
families. This is in contrast to integrable billiards (and to the 7r/2 internal angles in Fig. 1) 
where families of rays do not see the vertex and continue smoothly. 

We shall focus here on the periodic orbits of such systems for they form the skeleton 
on which generic classical motion is built. They are also the central objects of modern 
semiclassical theories || which provide a duality between the quantum spectrum and the 
classical length and stabilities of periodic orbits. In polygonal billiards, primitive periodic 
orbits can typically be classified under two categories depending on whether they suffer even 
or odd number of bounces at the boundary. In both cases however, the linearised flow is 
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marginally unstable since the Jacobian matrix, J p , connecting the transverse components 
(u±(t + T p ) = J p u±(t) where u± = (q±,p±) T and T p is the time period of the orbit) has unit 
eigenvalues. However, when the number of reflections, n p , at the boundary is even, the orbit 
occurs in a 1-parameter family while it is isolated when n p is odd. This follows from the fact 
that the left (right) neighbourhood of an orbit becomes the right (left) neighbourhood on 
reflection so that an initial neighbourhood can never overlap with itself after odd number of 
reflections ||. 

While isolated orbits are important and need to be incorporated in any complete periodic 
orbit theory, a family of identical orbits has greater weight and such families proliferate 
faster than isolated orbits . Besides, in a number of cases including the L-shaped billiard 
of Fig. 1, isolated periodic orbits are altogether absent since the initial and final momenta 
coincide only if the orbit undergoes even number of reflections. For this reason, we shall 
consider families of periodic orbits in this article (see Fig. (2) for an example). 

Not much is however known about the manner in which they are organised and the few 
mathematical results that exist |J concern the asymptotic properties of their proliferation 
rate. For a sub-class of rational polygons where the vertices and edges lie on an integrable 
lattice (the so called almost-integrable systems ||), these asymptotic results are exact. It is 
known for example that the number of periodic orbit families, N(l), increases quadratically 
with length, I, as I — > oo. For general rational polygons, rigorous results show that c\l 2 < 
< C2/ 2 for sufficiently large values of I |||7] while in case of a regular polygon P n with 
n sides, it is known that N(l) ~ c n l 2 / A || where c n is a number theoretic constant and A 
denotes the area of P n . Very little is known however about other aspects such as the sum 
rules obeyed by periodic orbit families in contrast to the limiting cases of integrable and 
chaotic behaviour where these have been well studied. Besides, it is desirable to learn about 
the variation of the proliferation rate as a function of the genus for this should tell us about 
the transition to chaos in polygonal approximations of chaotic billiards [pj, |TC)|1 . 

We shall concern ourselves primarily with a basic sum rule obeyed by periodic orbits 
arising from the conservation of probability. This leads to the proliferation law, N(l) = 
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Trbol 2 1 {a(l)) where 60 is a constant. The quantity (a(/)} is the average area occupied by all 
families of periodic orbits with length less than I and is not a constant unlike integrable 
billiards. We provide here the some numerical results on how (a(l)) changes with the length 
of the orbits and the genus of the invariant surface. While this does not allow us to make 
quantitative predictions, the qualitative behaviour sheds new light on the proliferation law 
for short orbits and its variation with the genus, g of the invariant surface. 

Finally, we shall also study correlations in the length spectrum of periodic orbits. The 
numerical results provided here are the first of its kind for generic pseudo-integrable billiards 
and corroborate theoretical predictions provided earlier ||11|| . 

The organisation of this paper is as follows. In section ||, we provide the basic sum 
rule and the proliferation law in pseudo-integrable systems. In section |HI| , we discuss al- 
gorithms to determine periodic orbit families in generic situations. This is followed by a 
numerical demonstration of the results obtained in section |I| as well an exploration of the 
area law. Correlations are discussed in sections [V] together with numerical results. Our 
main conclusions are summarised in section [VI]. 



II. THE BASIC SUM RULE 

The manner in which periodic orbits organise themselves in closed systems is strongly 
linked to the existence of sum rules arising from conservation laws. For example, the fact 



that a particle never escapes implies that for chaotic systems [fLlgJlS 



<r \ " T p 5{t rT p ) _ 

^^|det(l-J/)| ; 1 ] 

where the summation over p refers to all primitive periodic orbits, T p is time period, J p is the 
stability matrix evaluated on the orbit and the symbol (.) denotes the average value of the 
expression on the left. Since the periodic orbits are unstable and isolated, |det(l — J p r )| — 
e A P rT p, where X p is the Lyapunov exponent of the orbit. The exponential proliferation of 
orbits is thus implicit in eq. ([I]). 
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A transparent derivation of Eq. (JJ) follows from the classical evolution operator 

L*o0(x)= J dy 5(x - f(y)) 0(y) (2) 

where L l governs the evolution of densities </>(x), x = (q,p) and f* refers to the flow in 
the full phase space. We denote by A n (t) the eigenvalue corresponding to an eigenfunction 
n (x) such that L*o0 n (x) = A n (t)0 n (x). The semi-group property, L* 1 o L t2 = L tl+t ' 2 , for 
continuous time implies that the eigenvalues {A n (t)} are of the form {e An *}. Further, for 
hamiltonian flows, Eq. (§) implies that there exists a unit eigenvalue corresponding to a 
uniform density so that Ao = 0. For strongly hyperbolic systems, X n = —a n + i/3 n ,n > 1 
with a negative real part implying that 

Tr V = 1 + exp{-a n t + if3 n t} (3) 

n 

Eq. (HI) is thus a restatement of Eq. (|3|) with the trace expressed in terms of periodic orbit 
stabilities and time periods. 

For polygonal billiards, appropriate modifications are necessary to take account of the 
fact that the flow is restricted to an invariant surface that is two dimensional. Further, since 
classical considerations do not always yield the spectrum {\ n (t)}, we shall take resort to the 
semiclassical trace formula which involves periodic orbit sums similar to the kind that we 
shall encounter in the classical case. 

Before considering the more general case of pseudo-integrable billiards, we first introduce 
the appropriate classical evolution operator for integrable systems. This is easily defined as 

£*o0(0i, e 2 ) = J de' x de' 2 5(9, - o?)6{9 2 - o'i) 0' 2 ) (4) 

where 9\ and 8 2 are the angular coordinates on the torus and evolve in time as 9\ = 
co>j(ii, / 2 )£ + 9i with oji = dH(Ii,I 2 )/dIi and Zj = p.dq. Here Ti,i = 1,2 refer to 

the two irreducible circuits on the torus and p is the momentum conjugate to the coordi- 
nate q. 
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It is easy to see that the eigenfunctions {4> n {di,0 2 )} are such that <fin{0\, 0?) = 
A n (£)0 n (#i, 9 2 ) where A n (t) = e iant . On demanding that n (0i,02) be a single valued func- 
tion of (9 1 ,9 2 ), it follows that (f)n(di,9 2 ) = e % ( niei+n2e ^ where n = (ni,n 2 ) is a point on the 
integer lattice. Thus the eigenvalue, A n (t) = exp{it(nicji + n 2 u 2 )}. 

To illustrate this, we consider a rectangular billiard where the hamiltonian expressed in 
terms of the actions, ii, I 2 is H(Ii, I 2 ) = n 2 (If/L 2 + I 2 / 'L 2 ) where Li, L 2 are the lengths of 
the two sides. With I\ = \[EL\ cos(</?)/7r and I 2 = VEL 2 sin(</?)/7r, it is easy to see that at 
a given energy, i?, each torus is parametrised by a particular value of ip. Thus 

and the spectrum is continuous. The trace thus involves an integration over if as well as a 
sum over all ni, n 2 : 

Tr L ' = E rf V e l ^ sin( ^ n) = 2tt J] M\[e~J) (6) 

n J-K-iin. n 

where J is a Bessel function, / = 2t\/~E, tan(yU n ) = riiL 2 /(n 2 Li) and E n = it 2 (n\/L\ + 
n|/L|). On separating out n = (0,0) from the rest and restricting the summation to the 
first quadrant of the integer lattice, it follows that 

Tr L* = 27r + 27r7V^ J {JeJ). (7) 



where N = A. Note that the first term on the right merely states the fact that there exists a 
unit eigenvalue on every torus labelled by (p. Also, though we have not invoked semiclassics 
at any stage, the spectrum {E n } corresponds to the Neumann spectrum of the billiard 
considered. We shall subsequently show that this is true in general and for now it remains 
to express the trace of the evolution operator in terms of periodic orbits. 
The trace of L 1 in the integrable case can be expressed as 

Tr L l = Jd^J dQ x dQ 2 5{9 1 - 0*)<S(0 2 - 0*) (8) 

It follows immediately that the only orbits that contribute are the ones that are periodic. 
For a rectangular billiard, the integrals can be evaluated quite easily and yields 
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Tn« = 4EEr^«(i-W ( 9 ) 

Nl N2 lN lt N 2 

where {Ni,N 2 } are the winding numbers on the torus and label periodic orbits of length 
Ini,n 2 - Using Eqns. (|7|) and (|9|), it follows that 

4 E E - l NuN2 ) = 2vr + 2ixN E M^fKl) (10) 

Ni N2 N i> N 2 n 

Thus, the dominant non-oscillatory contribution to the trace that survives averaging is 2tt 
and this gives rise to the analogue of Eq. (JI|) 

4(EE^-W> = 2» (11) 

Ni N 2 l Ni,N 2 

The proliferation rate for the rectangular billiard thus follows from these considerations. 

These ideas can be generalised for polygonal billiards that are pseudo-integrable even 
though the structure of the invariant surface no longer allows one to use action and angle 
variables. 

For both integrable and pseudo-integrable polygonal billiards, the dynamics in phase 
space can be viewed in a singly connected region by executing 2g cuts in the invariant 
surface and identifying edges appropriately. At a given energy, the motion is parametrised 
by the angle, ip, that a trajectory makes with respect to one of the edges. As a trivial 
example, consider the rectangular billiard. The singly connected region is a larger rectangle 
consisting of four copies corresponding to the four directions that a trajectory can have and 
these can be glued appropriately to form a torus |TJ|]. As a non-trivial example, consider 



the L-shaped billiard of Fig. (1) which is pseudo-integrable with its invariant surface having, 
g = 2. Alternately, the surface can be represented by a singly connected region in the plane 
(see Fig. (3)) and consists of four copies corresponding to the four possible directions an 
orbit can have and these are glued appropriately. A trajectory in phase space thus consists 
of parallel segments at an angle (p measured for example with respect to one of the sides. 
It will be useful to note at this point that the same trajectory can also be represented by 
parallel segments at angles ir — ip, it + ip and 2ix — ip. In general, the number of directions 



for representing a trajectory equals the number of copies, N, that constitute the invariant 
surface. 

The classical propagator on an invariant surface parametrised by tp is thus 

L^)o0(q) = / dq' 5(q - q'^)) 0(q') (12) 

where q refers to the position in the singly connected region and q'*(<£>) is the time evolution 
parametrised by <p as described above. 

The trace of £*(y?) takes into account all possible invariant surfaces that exist and hence 
involves an additional integration over ip. Thus 

Tr 1} = J dtp J dq 5(q - q*(y?)) (13) 

Clearly the only orbits that contribute are the ones that are periodic. Further, the q integra- 
tions are simpler if we transform to a local coordinate system with one component parallel 
to the trajectory and the other perpendicular. Thus 5\\(q\\ — gjj) = \5(t — rT p ) where v is 
the velocity, T p is the period of the orbit and r is the repetition number. Similarly, for an 
orbit of period T p parametrised by the angle p p , 5±(q± — q±) — 8(ip — ( Pp)/\dq±/d<p\ (p=(fip 
where \dq±/d(p\ ip= = rl p for marginally unstable billiards. Putting these results together 
and noting that each periodic orbit occurs in general at N different values of p, we finally 
have 

Tr L* = £ A n (t) = iV £ E - rip) (14) 

n p r=l r V 

where I = tv and the summation over p refers to all primitive periodic orbit families with 
length l p and occupying an area a p . Note that Eq. (|) is a special case of Eq. (|14"D which 
holds for both integrable and non-integrable polygonal billiards. 

It is possible to re-express the periodic orbit sum in Eq. ( |14"D starting with the appropriate 



quantum trace formula [15 



1 oo 

£ 5(E - E n ) = d av (E) + -= £ ]T -j^= cos(krl p - - - vrrn p ) (15) 



^tt 3 p r=1 Jkrlp 4 



Here d av (E) refers to the average density of quantal eigenstates, k = y/E, l p is the length 
of a primitive periodic orbit family. The phase irrn p is set to zero while considering the 
Neumann spectrum while in the Dirichlet case, n p equals the number of bounces that the 
primitive orbit suffers at the boundary. For convenience, we have chosen h — 1 and the 
mass m = 1/2. Starting with the function 



I roo 

g(l) = £ /( sjE n l)e-^ = / dE f(VEl)e^ E £ 5{E - E n ) (16) 



where f(x) = y'^cos(x — 7r/4) and < e < E , it is possible to show using Eq. ( [To] ) that 
for polygonal billiards [17 



E E -£*(' - ^) = 27r& o + 2vr 2 /(\W) (17) 

p r=l ri P n 



for (3 — » + . In the above, 

a P (-l)^ 



&o = EE j ' f dEf(VEl)f(VErl p ) (18) 

„ r 4t Jo 



p r 



and is a constant []IT , |I8" |. It follows from Eqns. ( p~T| ) and (|H| ) that 

Tr L* = 53 E ~ = 27riY6 ° + 2ttJVE/(V^0- (19) 

p r =l r ^P n 

where {-E n } are the Neumann eigenvalues of the system. As in case of the rectangular 
billiard, the oscillatory contributions wash out on averaging so that 

(EE^-^)) = H (2°) 

This is the central result of this section and forms the basic sum rule obeyed by periodic 
orbit families. 

A few remarks about this derivation and the magnitude of bo are however in order. 
Eq. flT3| ) is exact for the L-shaped billiard and all other boundaries which preclude the 
existence of isolated periodic orbits. However, even for these shapes, the semiclassical trace 
formula is only an approximation to the exact density whenever the billiard in question 
is pseudo-integrable. In this sense, the sum rule in Eq. fl20|) is not expected to be exact. 
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However, we believe the existence of higher order corrections (such as diffraction) affects the 
magnitude of the constant bo while preserving the constant nature of the periodic orbit sum 
on the left. 

In the integrable case it is easy to show from other considerations that 60 = where 
N is the number of sheets that constitute the invariant surface 111 711. This also follows from 



Eqns. (|H]) and (|19D since boN = 1. For pseudo-integrable billiards, each invariant surface 
parametrised by ip has an eigenvalue, Ao(<^) = 1. Thus the non-oscillatory part of the trace 
should equal / Ao(y?) dip and to a first approximation this yields 2ir implying that 60 — 1/-/V. 
However each singular vertex connects two distinct points at any angle ip and hence the 
integration over ip is non-trivial. We can therefore state that 60 is approximately 1/N in the 
pseudo-integrable case while this is exact in the integrable case. We shall show numerically 
that the magnitude of deviations (from 1/N) in the pseudo-integrable case depends on the 
existence of periodic orbit pairs at the singular vertex. First however, we briefly describe 
the algorithms used to determine periodic orbit families. 

III. ALGORITHMS FOR DETERMINING PERIODIC ORBITS 

Periodic orbits in polygonal billiards are generally hard to classify. Unlike integrable 
billiards, they cannot be described by two integers which count the number of windings 
around the two irreducible circuits on the torus though in exceptional cases this can indeed 
be done. However, since the invariant surface has a well-defined genus, it is expected that a 
set of integers N = {Ni, N 2 , . . . , N 2g } obeying the relationship 

^i = -^r L (21) 

can be used to label periodic orbits. Here uji refers to the frequency corresponding to each 
irreducible circuit Tj and depends on the energy, E and the angle ip that labels each invariant 
surface. Note however that not all points on this multi-dimensional integer lattice are allowed 
since there are constraints and this method of labelling orbits becomes cumbersome for 
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surfaces of higher genus. Nevertheless, we illustrate the idea here for the L-shaped billiard 
of Fig. (1). 

Let the length of the two bouncing ball orbits in the X-direction be L\ and L 2 respectively 
and their lengths in the Y-direction be L3 and L4. These define the irreducible circuits for 
the L-shaped billiard. Thus : 

v cos(f) 2nNi 



2Li 




v cos((p) 


2irN 2 


2L 2 


T N 


vsm(<p) 


27r7V 3 


2L 3 


T N 


vsin(<f) 


2nN 4 


2L A 





(22) 



This implies that the angle (p at which a periodic orbit can exist is such that 



N 3 L 3 + N t L t 

Ua(V> = N lLl + N2L2 (23) 

Eqns. (|23|) and (^) merely express the fact that any periodic orbit should have integer 
number of windings around the irreducible circuits. Thus, the total displacement along the 
X-direction should be 2{N\L\ + N2L2) while the total displacement in the Y-direction should 
be 2(N 3 L 3 + N i L i ) where iVj are integers. As mentioned before however, not all realizations 
of {Ni} correspond to real periodic orbits and the final step consists in checking numerically 
whether a periodic orbit at the angle ip (given by Eq. (P5[)) exists. Note that one member 
of each family necessarily resides at one of the singular vertices and hence it is sufficient to 
verify the existence of this orbit. 

This method works equally well for other billiards with steps (see Fig. (1)) and the 
number of integers necessary to describe orbits increases with the number of steps. An 
alternate method which exploits the fact that periodic orbits occur in families is often useful 
when the irreducible circuits are not obvious and this is described below. 

Note that a non-periodic orbit originating from the same point q (e.g. the singular 
vertex) as a periodic orbit but with a momentum slightly different (from the periodic orbit) 
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suffers a net transverse deviation q± = (— l) n ^ sm(ip—(pp)l,p. Here l v is the distance traversed 
by a non-periodic orbit at an angle <p after n v reflections from the boundary and (p p is the 
angle at which a periodic orbit exists. This provides a correction to the initial angle and a 
few iterations are normally sufficient to converge on a periodic orbit with good accuracy. In 
order to obtain all periodic orbits, it is necessary to shoot trajectories from every singular 
vertex since one member of each family resides at one of these vertices. 

Apart from the length of a periodic orbit, it is also important to compute the area 
occupied by the family. This can be achieved by shooting a single periodic trajectory which 
resides at a singular vertex and by noting that this orbit lies on the edge of a family. 
Thus if the rest of the family lies on the left neighbourhood initially, it is necessary to 
determine the perpendicular distance from a singular vertex to this trajectory every time the 
initial neighbourhood lies towards this singular vertex. The shortest of these perpendicular 
distances gives the transverse extent of the family and the area can thus be computed. 

These algorithms have been used to generate the lengths {l p } and areas {a p } of primitive 
orbits in the L-shaped billiard as well as the two and three step billiards. We present our 
numerical results in the following sections. 



IV. NUMERICAL RESULTS : THE SUM RULE AND AREA LAW 

We present here our numerical results on the 1-step (L-shaped), 2-step and 3-step bil- 
liards. Their invariant surfaces have genus 2, 3 and 4 respectively and the quantities we 
shall study are the sum rule derived in section [II] and the variation of (a(l)) with I and the 
genus, g of the invariant surface. 

The sum rule we wish to study can be re-expressed as 

SW^EE £=W (24) 

and this is plotted in Fig. (4) for four different 1-step (L-shaped) billiards. Notice first that 
in each case the behaviour is linear as predicted by Eq. (E4f). Besides, in three of the four 
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cases, the slopes are quite close (6 — 0.27) and these correspond to non-degenerate L-shaped 
billiards with sides that are unequal and irrationally related. The one with a substantially 
larger slope (bo ~ 0.32) is a degenerate case of a square with a quarter removed and for 
this system there exist a substantial number of periodic orbit pairs at the same angle on 
the two adjacent edges at the singular vertex. The differences between the degenerate and 
non-degenerate case also persists in other quantities as we shall shortly see. 

We next plot S(l) for non-degenerate examples of a 1,2 and 3-step billiards in Fig. (5). 
The number of orbits considered are far less due to increased computational effort though 
the linear behaviour is obvious in all three cases. The slopes are again close to 0.25 and vary 
from b ~ 0.24 to b ~ 0.28. 

Thus, periodic orbits obey a basic sum rule given by Eq. fl24|) where bo is a constant. 
Further, the magnitude of bo is close to 0.25 in all cases and the deviations from this value 
are larger in the degenerate case where periodic orbit pairs exist at the singular vertex. 

The sum rule yields the proliferation rate of periodic orbits 

m) = (25) 

where (a(l)) is the average area occupied by periodic orbit families with length less than /. 
In case of integrable polygons (a(l)) is a constant and equals AN where A is the area of the 
billiard while b = l/N. For pseudo- integrable cases, the areas a p occupied by individual 
periodic orbit families generically occupy a wide spectrum bounded above by NA. As the 
length of a family increases, encounters with singular vertices are more frequent so that the 
transverse extent of family decreases. Quantitative predictions about the behaviour of (a(l)) 
are however difficult and we provide here some numerical results. 

Fig. (6) shows a plot of [a(l))/NA for three different L-shaped billiards one of which 
is the degenerate case presented in Fig. (4). The average area increases initially before 
saturating for large I in all cases. The normalised saturation value seems to be independent 
of the ratio of the sides as long as they are irrationally related (we have observed this for 
other cases not presented here) but is very different for the degenerate example of a square 
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with a quarter removed. 

A comparison of non-degenerate 1,2 and 3-step billiards is shown in Fig. (7). Recall that 
the invariant surfaces for these have genus respectively equal to 2, 3 and 4. The saturation 
observed for the 1-step (L-shaped) case seems to be a common feature and interestingly the 
normalised average saturation value decreases. This is however expected since the number 
of singular vertices increases with the genus thereby reducing the average transverse extent 
of orbit families at any given length. 

These observations allow us to conclude that for short lengths, the proliferation law is 
sub-quadratic. Asymptotically however, N(l) ~ I 2 as in integrable billiards. Further, the 
asymptotic proliferation rate for billiards with the same area A increases with the genus due 
to a decrease in the asymptotic normalised saturation value of (a(l)). 

These numerical results provide a qualitative picture of the area law and show that peri- 
odic orbits in polygonal billiards are organised such that they obey a sum rule. Quantitative 
predictions would require a more extensive numerical study such that empirical laws for the 
saturation value and its variation with the genus can be arrived at. 

V. CORRELATIONS IN THE LENGTH SPECTRUM 

Our numerical explorations so far have been focussed on the average properties of the 
length spectrum. We shall now attempt to understand the nature of the fluctuations and 
characterise their statistical properties. 

Fluctuations in the length spectrum are generally difficult to study from purely classi- 
cal considerations. There are notable exceptions however. Fluctuations in the integrable 
case can be studied using the Poisson summation formula since the lengths of orbits are 
expressed in terms of integers {A 1; N 2 } ||19| . The other extreme is the motion on surfaces 
of constant negative curvature where the Selberg trace formula provides an exact dual rela- 
tionship between the classical lengths and the eigenvalues of the Laplace-Beltrami operator 
PH] . For other systems, a possible way of studying fluctuations in the length spectrum lies 
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in inverting the semiclassical quantum trace formula. For pseudo-integrable billiards, this 
has been achieved in section [TT] and the integrable density of lengths, N(l) can be expressed 
as 

m = + -rjjSS 12 f dl ' l'f(\fEj) (26) 

(a{l)) (a{l)}^Jo v 

Statistical properties of the fluctuations can thus be studied using techniques introduced for 
the fluctuations in the quantum energy spectrum [ 2l|,|TT . 



The correlations commonly studied are the nearest neighbour spacings distribution, P(s) 
and a two-point correlation referred to as the spectral rigidity, As(L). The rigidity measures 
the average mean square deviation of the staircase function N(l) from the best fitting straight 
line over L mean level spacings. For a normalised (unit mean spacing) Poisson spectrum, 
P{s) = e~ s while A 3 (L) = L/15. These features are commonly observed in the quantum 
energy spectrum of integrable systems as well as in their length spectrum. For chaotic 
billiards, the quantum energy spectrum generically produces non-Poisson statistics while 



the length spectrum correlations are Poisson for long orbits at least over short ranges p2 
There are however deviations that can be observed in A^(L) for short orbits or over longer 
ranges in the spectrum J2"2|. With this background, we now present our results for pseudo- 



integrable billiards. 

Figs. (8) and (9) show plots of P{s) and A 3 (L) for a non-degenerate L-shaped billiard. 
A total of 3000 lengths have been considered after excluding the shortest 3000 orbits. The 
correlations are clearly Poissonian as in case of integrable or chaotic billiards. 

For the 3-step billiard where fewer lengths are available (about 1250), we have carried 
out a similar study and the results are shown in Figs. (10) and (11). Deviations from the 
Possion behaviour can now be seen especially in the spectral rigidity. By considering shorter 
lengths in the L-shaped billiard, similar deviations were observed. 

The statistical properties of fluctuations in the length spectrum of pseudo-integrable 
systems are thus similar to those of chaotic billiards. For long orbits the correlations are 
Poisson while deviations exist for shorter orbits. 
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VI. DISCUSSIONS AND CONCLUSIONS 



Quantum billiards are experimentally realizable in the form of microwave cavities since 



the wave equations are identical [23,23|. They are also interesting in their own right and 
have proved to be the testing grounds for several ideas in the field of Quantum Chaos. Of all 
possible shapes, polygonal billiards have perhaps been the least understood largely because 
very little was known about the organisation of periodic orbits. The results presented 
here have however been used recently to obtain convergent semiclassical eigenvalues with 
arbitrary non-periodic trajectories |25| as well to demonstrate that periodic orbits provide 
excellent estimates of two-point correlations in the quantum energy spectrum |1| . Though it 



is beyond the scope of this article to review these recent developments on quantum polygonal 
billiards, we refer the interested reader to these articles and the references contained therein. 

In the previous sections we have explored the organisation of periodic orbits in rational 
polygonal billiards that are pseudo-integrable. Our main conclusions are as follows : 

• Orbit families obey the sum rule {J2 P J2^i a p 8{l — rl p )/{rl p )) = 27r6 thereby giving 
rise to the proliferation law N(l) = irbol 2 / (a(l)) for all rational polygons. 

• The quantity bo is approximately 1/N in generic generic rational billiards and deviations 
from this value are observed to be significant in degenerate situations. 

• (a(l)) increases initially before saturating to a value much smaller than the maximum 
allowed area NA. The asymptotic proliferation law is thus quadratic even for systems that 
are not almost-integrable and the density of periodic orbits lengths is far greater than an 
equivalent integrable systems having the same area. 

• The normalised average area (a(l))/NA decreases with the genus of the billiard while 
60 is approximately 0.25. Periodic orbits thus proliferate faster with an increase in genus. 

• The statistical properties of the fluctuations in the length spectrum are Poisson when 
the orbits considered are long. 
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FIGURES 



Y 



X 

FIG. 1. The 1, 2 and 3-step billiards are pseudo-integrable due to the 3ir/2 internal angles. The 
1-step billiard (bottom) is also referred to as the L-shaped billiard in the text 
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FIG. 2. An example of a periodic orbit in the 3-step billiard. Note that adjacent parallel rays 
also periodic so that a family of periodic orbits exists. Its extent is limited by the 3ir/2 vertices. 
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<?1 

FIG. 3. The singly connected region for an L-shaped billiard consists of four copies with edges 
appropriately identified. A trajectory originating near the 3ir/2 vertex in 1 is plotted in configu- 
ration space using bold lines and the corresponding unfolded trajectory is also shown. The latter 
consists of parallel segments and the trajectory can be parametrized by the angle ip that it makes 
for example with the q\ axis. 
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S(l) 




FIG. 5. The sum rule for a 1,2 and 3-step billiard. They are respectively the middle, bottom 
and top curves. 
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FIG. 6. The normalized average area, («(/)) = (a(l)) /N A as a function of the length / for 3 
examples of the 1-step billiard. The saturation value differs for the degenerate example mentioned 
earlier 
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FIG. 7. The normalized area for an example each of the 1-step (top), 2-step (middle) and 3-step 
(bottom) billiards. 
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FIG. 8. The nearest neighbour spacings distribution for a non-degenerate example of a 1-step 
billiard. Also plotted is the Poisson distribution. 
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